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Abstract. We calculate the structure of accretion disk around a spinning black hole for accretion 
rates M = 0.01 — 10M a s . The model is fully relativistic and treats accurately the disk micro- 
physics including neutrino emissivity, opacity, electron degeneracy, and nuclear composition. We 
find that the accretion flow always regulates itself to a mildly degenerate state with the proton- 
to-nucleon ratio Y e ^0.1 and becomes very neutron-rich. The disk has a well defined "ignition" 
radius where neutrino flux raises dramatically, cooling becomes efficient, and Y e suddenly drops. 
We also calculate other characteristic radii of the disk, including the v-opaque and v-trapping radii, 
and show their dependence on M. Accretion disks around fast-rotating black holes produce intense 
neutrino fluxes which may deposit enough energy above the disk to generate a GRB jet. 
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Most GRB models assume a short-lived accretion disk as a primary source of energy. Such 
disks have huge accretion rates M = 0.01 — 1OM s~ l and their cooling mechanism is neutrino 
emission flylS]- Special features of neutrino-cooled disks are: 

— The main cooling is due to reactions p + e~ — » n + v and n + e + — > p + v. 

— Electron degeneracy affects the disk structure. 

— The disk can be opaque to neutrinos which affects the cooling rate. 

— Nuclear composition changes with radius, consuming energy. 

— Neutrino cooling and electron degeneracy lead to very high neutron richness. 

The disk model may be formulated following the classical vertically-integrated approach of 
Shakura & Sunyaev and approximating the azimuthal motion in the disk, u^, as Keplerian 
rotation. A small radial velocity is superimposed on this rotation. It is approximately given by 
u r = aB(H /r) 2 ru^ where H is half-thickness of the disk (found from hydrostatic balance), B is 
a numerical factor determined by Kerr metric, and a = 0.1 — 0.01. 

The main equations of the accretion disk are as follows. 



Energy conservation : 



cu 



d(UH) U + Pd(pH) 



dr 



dr 



(1) 



Lepton conservation : 
Charge conservation : 

Chemical balance : 
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Here F + is the viscously dissipated energy per unit area, F is the neutrino energy flux, N v and 
Ny are the number fluxes of neutrinos and antineutrinos, p is the mass density, T is temperature, 
Y e is the proton-to-baryon ratio, jx e and jx v are the chemical potentials of electrons and neutrinos, 
n e - and n e + are number densities of electrons and positrons. The energy density U and pressure 
P include contributions of baryons, electrons, positrons, radiation, and neutrinos. They are 
calculated by integrating the corresponding distribution functions (Fermi-Dirac distributions for 
e* 1 and V, V if the disk is opaque to neutrinos). 

Note that we include the advection terms in the laws of energy and lepton conservation (right- 
hand side of eqs. 1 and 2). They describe the most important effects of advection. Other effects 
of heat storage in the disk — in particular, the dynamical effects of the radial pressure gradient 
on u r and — are relatively small and neglected here, i.e. the standard a-disk model is used to 
calculate vfi and u r . 

The chemical balance and jU v are used only in the v-opaque region; then the cooling rate is 
calculated as F = cU v / T v + cU v /ty. In the v-transparent region, F equals 2H(U e - p + U e + n ) 
where U e - p and U e + n are the rates of energy release by reactions e~ + p — ► n + v and e + + n — > 
p + V. The abundance of a-particles is determined by the nuclear statistical equilibrium. 

We have solved the set of disk equations numerically starting from an outer region where 
no neutrino cooling takes place and assuming that the infalling matter is initially made of a 
particles (heavier elements are decomposed into a-particles as the matter approaches the black 
hole). A schematic picture of the accretion disk and its characteristic radii are shown in Figure 1. 
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FIGURE 1. Schematic picture of the disk with characteristic radii indicated. 

The outer region r > 50r g (r g = 2GM/c 2 is the gravitational radius of the black hole) is dom- 
inated by a particles. Neutrino cooling is negligible in this region and the viscously dissipated 
heat is stored in the disk and advected inward. The destruction of a particles at r a ~ 50r g con- 
sumes ~ 7 MeV per baryon, which makes the disk somewhat thinner. At the "ignition" radius 
r ign (slightly smaller than r a ), kT reaches ~ MeV and neutrino cooling becomes significant, 
further reducing disk scale-height H/r. With decreasing radius the disk becomes opaque to neu- 
trinos (first to v and then to v). In the innermost region of disks with M > 2M s _1 the neutrinos 
are trapped and advected toward the black hole. Figure 2 shows the found characteristic radii as 
functions of M for disks with a = 0.1 around black holes with spin parameters a = and 0.95. 

We find two general properties of the disk: 
1. — It regulates itself to a mildly degenerate state with \i e = 1 ~ 3kT. The reason of this 
regulation is the negative feedback of degeneracy and neutronization on the cooling rate: higher 
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FIGURE 2. Characteristic radii in the accretion disk as a function of M for a black hole of mass 
M = 3Mq and spin parameter a = 0.95 (left) and (right). Viscosity parameter a = 0.1 is assumed. 
Neutrino cooling does not ignite when M is below 0.02Mqs~' for Kerr and 0.1M Q ,s for Schwarzschild 
black holes. Neutrino trapping takes place for Kerr black holes when M > 2M Q s _1 . The disk extends 
down to the marginally stable orbit, which is 3r ? for a = and as r g for a — 0.95. 

degeneracy }X e /kT —* fewer positrons (n + /n_ ~ e~^' kT ) —> lower equilibrium rate of neutrino 
emission — > lower cooling rate — ► higher temperature — > lower degeneracy. 
2. — In a broad range of relevant accretion rates the disk is found to reach Y e ~ 0.1. It 
corresponds to a very high neutron-to-proton ratio ~ 10. 

Figure 3 shows examples of the disk structure around a Kerr black hole with a = 0.95 with 
accretion rate O.2M s _1 , for three different values of viscosity parameter a. The transitions 
at the characteristic radii manifest themselves in the non-monotonic behavior of the degeneracy 
parameter r] = jj, e /kT, F /F + ', Y e , and// Jr. In particular the sharp maximum of F~ /F + appears 
at the ignition radius. The disk is quite thin at r < rj gn and cooled locally, F « F + . 

The high neutron richness has important implications for the global picture of GRB explosion. 
When the neutron-rich material is ejected in a relativistic jet and develops a large Lorentz factor, 
the neutrons decay at distance ~ 10 16 cm and affect the observed explosion. 

Disks around rapidly rotating black holes produce a high neutrino flux. It deposits energy 
above the disk (via neutrino annihilation v + V — » e + + e~) and can drive an outflow with power 
L > 10 51 erg/s. A detailed study of vv annihilation around a Kerr black hole is in preparation. 
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FIGURE 3. Numerical solution for accretion disk with M = O.2M ,s around a black hole with mass 
M = 3Mq and spin parameter a = 0.95. Three models are shown with viscosity parameter a = 0.1, 0.03 
and 0.01. The degeneracy parameter r\ is defined as jj, e /kT, Xf ree is the mass fraction of free nucleons. 



